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lar  Materials  and  Lubrication  Phenomena,”  Task  No.  734203,  “Fundamental  Principles 
Determining  the  Behavior  of  Macromolecules.”  It  was  administered  under  the  direction 
of  the  AF  Materials  Laboratory,  Research  and  Technology  Division,  Lt.  Martin  M.  Tessler, 
project  engineer. 
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ABSTRACT 


The  random  thermal  degradation  of  four-  and  six-membered  ring  ladder  pol5miers 
were  investigated  and  results  compared  to  a  single  chain  polymer  undergoing  degradation 
under  the  identical  conditions.  The  change  in  molecular  weight  versus  time  was  plotted 
and  significant  differences  were  found  in  the  shapes  of  the  curves. 

Similar  studies  were  conducted  with  imperfect  ladder  polymers  which  had  single  bonds 
along  the  backbone  of  the  polymer. 


This  technical  documentary  report  has  been  reviewed  and  is  approved. 


WILLIAM  E.  GIBBS 
Chief,  Polymer  Branch 
Nonmetallic  Materials  Division 
AF  Materials  Laboratory 
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INTRODUCTION 


Interest  in  the  development  of  thermally  stable  materials  has  produced  much  research 
into  new  and  novel  types  of  polymers.  One  area  of  investigation  is  the  preparation  of  a 
ladder  or  double-strand  polymer.  The  first  successful  synthesis  of  a  ladder  polymer  was 
reported  by  Brown  and  co-workers  (Reference  1).  They  prepared  the  double-chain 
polyphenylsilesquioxane  in  which  the  polymer  backbone  is  completely  inorganic.  Attempts 
to  synthesize  a  ladder  polymer  with  an  organic  backbone  from  vinyl  isocyanate  (Refer¬ 
ence  2)  poly-3,4-isoprene  (Reference  3)  and  conjugated  dienes  (Reference  4)  have  recently 
been  reported.  These  polymers  have  segments  of  fused  rings,  but  they  are  not  completely 
fused  into  the  desired  ladder  structure. 

This  report  compares  the  thermal  stability  of  a  ladder  polymer  undergoing  random 
degradation  with  that  of  a  single  chain  polymer.  A  digital  computer  is  used  to  set  up  a 
statistical  (Monte  Carlo)  (Reference  5)  model  of  the  degrading  system.  The  following 
assumptions  are  made  in  defining  the  degrading  system: 

1.  The  polymer  sample  is  initially  homogeneous;  that  is,  only  chains  of  one  length  are 
present. 

2.  All  bonds  in  the  polymer  chain  are  of  equal  strength  and  equal  accessibility,  re¬ 
gardless  of  their  positions  in  the  molecule  and  the  length  of  the  chain. 

3.  The  rate  of  breaking  of  bonds  is  proportional  to  the  number  of  bonds  present  in  the 
degrading  system, 

4.  The  system  is  closed;  that  is,  no  fragments  can  leave  the  system. 

The  ratio  of  the  number  average  molecular  weight  (at  time  kt)  to  the  initial  number 
average  molecular  weight  is  obtained  directly  from  the  computer,  where  k  is  the  pro¬ 
portionality  constant  between  the  rate  of  breaking  bonds  and  the  number  of  bonds  present 
in  the  system  at  time  t. 


MODEL  LADDER  STRUCTURES 


Two  tjqjes  of  ladder  polymers  were  studied.  Type  I  is  a  fused  four-membered  ring 
and  Type  II  is  a  fused  six-membered  ring.  All  the  conclusions  reached  for  a  four- 
membered  ring  are  equally  applicable  to  an  eight- membered  ring. 


Manuscript  released  by  the  author  April  1964  for  publication  as  an  RTD  Technical 
Documentary  Report. 
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Initially,  a  simplified  scheme  of  degradation  was  considered.  In  this  simple  case,  a 
broken  bond  can  break  a  molecule  only  if  the  lx)nd  opposite  it  is  broken.  The  fact  that 
a  broken  crosslink  will  increase  the  probability  of  a  molecule  breaking  is  ignored.  This 
type  of  degradation  is  shown  in  Figure  1.  The  complex  case  of  degradation  is  shown  in 
Figure  2  where  breaks  in  the  crosslinks  will  result  in  increased  molecule  breaks.  This 
is  a  much  more  realistic  picture  of  random  degradation  than  the  simple  case. 


TYPE  I  TYPE  n 


Figure  1.  Different  Ways  a  Molecule  Can  Break  -  Simple  Case 


TYPE  H 

Figure  2.  lEfferent  Ways  a  Molecule  Can  Break  -  Complex  Case 


MONTE  CARLO  MODEL 


To  construct  a  Monte  Carlo  model  for  the  degradation  reaction,  a  portion  of  the  com¬ 
puter  storage  is  set  aside  to  represent  the  polymer  molecules.  Each  storage  location 
represents  one  bond.  If  a  bond  is  broken,  the  number  one  is  stored  in  its  location  and  if 
it  is  not  broken,  a  zero  is  stored  there.  Initially,  all  the  storage  locations  are  set  to  zero. 
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A  random  number  is  generated  which  has  the  limits  of  O  and  Bq-1,  where  Bq  is  the 

total  number  of  bonds  initially  present  in  the  system.  Each  number  then  represents  a 
particular  bond  in  a  particular  molecule  and  the  number  of  molecules  and  the  number  of 
bonds  per  molecule  can  be  varied  at  will  with  the  limiting  factor  being  the  size  of  the 
computer  storage.  In  the  present  study,  a  system  of  100  molecules  with  199  bonds  per 
molecule  for  ladder  polymer  Type  I  and  100  molecules  with  196  bonds  per  molecule  for 
ladder  polymer  Type  II  were  chosen.  Simultaneously,  a  single  chain  polymer  with  the 
same  number  of  bonds  per  molecule  was  degraded  in  the  computer  for  comparison 
purposes. 

When  a  random  number  is  generated,  the  bond  which  it  represents  is  checked  to  see 
if  it  is  broken  or  not.  If  it  is  already  broken,  a  duplicate  is  recorded  and  a  new  random 
number  generated.  If  it  is  not  broken,  the  bond  is  set  equal  to  one  and  the  opposite  bond 
checked.  If  the  opposite  bond  is  broken,  a  broken  molecule  is  recorded  and  a  new  random 
number  is  then  generated.  If  the  opposite  bond  is  not  broken,  no  further  work  is  done  in 
the  simple  case  and  a  new  random  number  is  generated.  In  the  complex  case,  the  com¬ 
puter  checks  all  the  crosslinks  to  the  left  and  right  until  it  either  comes  to  an  unbroken 
crosslink  or  the  last  crosslink  at  the  end  of  the  molecule.  If  the  last  crosslink  is  broken, 
a  broken  molecule  is  recorded  and  a  new  random  number  generated.  If  the  computer  finds 
an  unbroken  crosslink,  it  then  scans  all  of  the  side  chain  bonds  on  both  sides  of  the  ladder 
pol3nmer  from  the  broken  bond  produced  by  the  random  number  generation  to  the  last  side 
chain  bond  before  the  unbroken  crosslink.  If  any  of  the  side  chain  bonds  are  broken,  a 
broken  molecule  is  recorded  and  a  new  random  number  is  then  generated.  If  none  of  the 
side  chain  bonds  are  broken,  a  new  random  number  is  generated.  The  computer  keeps 
track  at  all  times  of  the  number  of  duplicates  recorded.  The  number  of  bonds  which  are 
broken  is  equal  to  the  number  of  random  numbers  generated  minus  the  number  of 
duplicates. 

The  number  of  molecules  present  in  the  system  (at  time  kt)  is  equal  to  the  number  of 
molecules  initially  present  plus  the  number  of  broken  molecules  since  each  broken 
molecule  produces  one  more  fragment.  The  time  factor  (kt)  is  equal  to  In  (Bq/B)  where 

Bq  is  the  number  of  bonds  present  initially  and  B  is  the  number  of  bonds  present  at  time 

kt.  The  ratio  of  the  number  average  molecular  weight  at  time  kt,  to  the  initial  number 
average  molecular  weight  is  equal  to  100/(100  +  number  broken  molecules).  All  of  these 
quantities  are  readily  calculated  during  the  course  of  the  degradation  and  the  computer 
prints  out  the  desired  data  at  any  convenient  interval  of  kt. 

The  subprogram  by  which  the  computer  generates  a  random  number  makes  use  of  an 
equation  in  which  multiplication  and  division  are  used,  the  output  being  the  remainder 
after  division.  Several  different  sets  of  random  numbers  were  used  to  check  the  accuracy 
of  the  method.  The  random  number  generator  appears  to  be  very  accurate  and  reliable. 


RESULTS  FOR  PERFECT  LADDER  POLYMERS 


The  curve  for  the  ratio  of  the  number  average  molecular  weight  at  time  kt  to  the  initial 
number  average  molecular  weight  versus  kt  is  shown  for  ladder  polymer  Type  I  in 
Figure  3  and  for  ladder  polymer  Type  II  in  Figure  4.  Both  the  simple  case  and  the  com¬ 
plex  case  are  shown  as  well  as  the  analogous  single  chain  polymer.  The  proportionality 
constant  k  is  not  known  and  it  probably  will  have  different  values  for  different  polymers. 
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This  will  shift  the  curves,  but  the  sigaificant  fact  is  the  shapes  of  the  curves.  The  single 
chain  polymer  shows  a  very  sharp  drop  in  molecular  weight  as  soon  as  degradation  begins. 
The  ladder  polymer  has  an  induction  period  where  the  molecular  weight  hardly  changes 
and  then  it  drops  much  more  slowly  than  the  corresponding  sin^e  chain  polymer.  As 
degradation  proceeds,  both  curves  merge  together.  These  results  indicate  that  the  ladder 
polymers  should  have  increased  thermal  stability  over  single  chain  polymers,  but  the 
extent  of  this  stability  can  only  be  determined  by  studying  the  chemistry  of  these  compounds. 


IMPERFECT  LADDER  POLYMERS 


Ladder  polymers  are  of  interest  because  of  their  potential  stability  at  high  temperatures. 
The  double-chain  polyphenylsilsesquioxane,  synthesized  by  Brown  and  co-workers 
(Reference  1)  has  an  inorganic  ladder  structure  and  the  silicon-oxygen  bonds  rearrange 
at  elevated  temperatures.  The  reported  attempts  to  sjmthesize  ladder  polymers  with  an 
organic  ladder  structure  (References  2,  3  and  4)  have  all  used  the  same  general  approach. 

A  polymer  is  first  prepared  which  has  pendant  nonconjugated  1,6  unsaturation  along  the 
chain  of  the  molecule.  Overberger  and  co-workers  (Reference  2)  polymerized  vinyl 
isocyanate  through  the  vinyl  group  or  the  isocyanate  group  to  get  the  starting  polymer, 
Angelo  (Reference  3)  prepared  poly-3,4-isoprene  and  Gaylord  and  co-workers  (Refer¬ 
ence  4)  used  butadiene,  isoprene  and  chloroprene  to  obtain  the  desired  polymer.  Cycli- 
zation  of  the  pendant  nonconjugated  1,6  unsaturated  polymer  with  suitable  catalysts 
produces  a  ladder  structure. 

The  problems  in  preparing  a  completely  perfect  ladder  polymer  are  enormous. 

Imperfect  ladder  polymers  can  result  from  isolated  unreactive  groups,  from  alternative 
cyclizations  or  isomerizations,  from  crosslinking  and  from  chain  scission  reactions. 
Extremely  dilute  solutions  must  be  used  to  prevent  crosslinking.  Another  problem  is 
that  as  the  ladder  polymer  is  being  synthesized,  it  solubility  is  decreasing  rapidly  and 
it  may  precipitate  out  of  solution  before  ring  closure  is  complete.  If  the  polymer  with 
the  pendant  nonconjugated  1,6  unsaturation  is  not  strucmrally  perfect,  complications 
will  arise  in  the  cyclization  step.  For  example,  the  poly- 3,4 -isoprene  used  by  Angelo 
contained  up  to  10  percent  of  1,4  structural  units. 

The  difficulty  in  preparing  perfect  double-chain  polymers  raises  the  question  as  to 
how  the  thermal  stability  of  an  imperfect  ladder  polymer  compares  with  that  of  a  single 
chain  polymer.  A  partial  solution  to  this  problem  can  be  obtained  by  setting  up  in  the 
storage  of  a  digital  computer,  an  imperfect  ladder  polymer.  Each  storage  location  will 
correspond  to  a  different  bond  and  a  random  number  generator  will  be  used  to  locate  the 
bond  which  is  being  broken  during  degradation.  An  unbroken  bond  will  be  represented  by 
a  zero  and  a  broken  bond  by  a  one.  If  all  the  bonds  are  initially  set  equal  to  zero,  we 
will  have  a  perfect  ladder  polymer.  If  we  initially  set  a  bond  equal  to  one  at  prescribed 
intervals  along  the  chain  we  will  have  an  imperfect  ladder  polymer  with  the  “broken 
bonds”  corresponding  to  missing  bonds  in  an  imperfectly  synthesized  ladder  polymer. 

If  the  random  number  generator  should  produce  a  number  corresponding  to  one  of  these 
missing  bonds,  a  duplicate  will  be  recorded  and  a  new  random  number  will  be  generated 
without  any  further  calculations.  It  is  assumed  that  the  polymer  sample  is  homo¬ 
geneous,  ail  bonds  in  the  polymer  have  an  equal  probability  of  being  broken,  the  rate  of 
breaking  bonds  is  proportional  to  the  number  of  bonds  present  in  the  degrading  system 
and  no  fragments  can  leave  the  system. 
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The  Monte  Carlo  calculations  are  identical  to  those  described  previously. 

Four-  (Type  I)  and  six-  (T5^e  II)  membered  ring  ladder  polymers  were  studied  and 
both  the  simple  and  complex  cases  of  degradation  examined.  A  single  chain  polymer 
with  the  same  number  of  bonds  per  molecule  was  degraded  simultaneously  in  the  com¬ 
puter  for  comparison  purposes. 

Type  1  ladder  polymers  were  studied  with  single  bonds  first  at  every  sixth  and  then 
at  every  twelfth  bond  along  the  chain.  Type  II  ladder  polymers  were  studied  with  single 
bonds  first  at  every  twelfth  and  then  at  every  twenty-fourth  bond  along  the  chain. 


RESULTS  FOR  IMPERFECT  LADDER  POLYMERS 


The  introduction  of  single  bonds  along  the  chain  of  the  ladder  polymer  results  in  a 
sharp  decrease  in  molecular  weight,  as  degradation  occurs,  compared  to  a  perfect 
ladder  polymer.  The  stability  of  the  imperfect  ladder  polymer  is  still  much  greater  than 
that  of  an  analogous  single  chained  polymer.  The  shape  of  the  curves  when  the  ratio  of 
the  number  average  molecular  weight  at  time  kt  to  the  initial  number  average  molecular 
weight  is  plotted  against  kt  is  very  interesting.  With  a  small  number  of  single  bonds 
within  the  ladder  chain,  the  curves  still  have  an  induction  period  where  the  change  of 
molecular  weight  is  relatively  small.  It  then  drops  very  rapidly,  although  not  as  fast 
as  a  single  chain  polymer.  As  the  number  of  sin^e  bonds  are  increased,  the  induction 
period  disappears  and  the  shape  of  the  curve  is  the  same  as  that  of  the  single  chain 
polymer  although  the  molecular  weight  still  does  not  decrease  as  fast.  With  a  suitable 
value  of  k,  it  is  possible  that  the  curves  will  become  superimposable  and  that  the  im¬ 
perfect  ladder  polymer  will  have  no  increased  advantage  in  thermal  stability  over  a 
single  chain  polymer.  The  answer  to  this  problem  must  await  the  synthesis  and  charac¬ 
terization  of  imperfect  ladder  polymers  and  a  study  of  their  chemical  behavior  when 
undergoing  thermal  degradation. 
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APPENDIX  I 
COMPUTER  PROGRAM 
SIMPLE  CASE.  TYPE  I 
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C  TESSLER  LAOOER  NJLYMER  DEGRADATION 
DIMENSION  LI 100,200) 

FREOUENCY  71 1 ,0,2 ) »8 1 9, 1 ,0 ) ,910 , 1 ,9 1 *4l 1 ,0 » 1 ) ,10 (9»1 ,0) , 13 (0,1 ,9) 
1,14I9»1*0) 

NBRAK»0 

NDUP»0 

NLNK=0 

DO  25  IX=1»30 
DO  24  1=1*200 
R=RAN01IY 1*19900, 

NR=R 

ML«NR/199  +1 
NB*199*ML  -  NR 
7  IF|NB-56)8,8,4 

8  IFrLIML*NB1-l)28,il*2« 

28  L(ML,NB)=1 

9  IFIl-LIML,NB+66)i24,  5*24 

4  1F(NB-132 110,10,14 

10  1F|LCML*NB»-1J17,11,17 
17  LCML,NB)»1 

13  IF<l-L(ML,NB-86n24,5*24 

14  IF(L(ML*NB)>1)12,11,12 

12  L(ML*N8)»1 

NLNK»NLNK+i 
GO  TO  24 

5  NBRAK=NBRAK  +1 
GO  TO  24 

11  NDUP=M0UP+1 

24  CONTINUE 
BRK=NBRAK 

ANMW=iOO, /< 100,+BRK ) 

X!  =  IX 
DUPN=N0UP 

T 1 ME-L06F 1 19900 , / f 19900 , -200 ,*X 1 +DUPN I ) 
ASMW=100,/»100,+200,*XI-DUPN) 

WRITE  OUTPUT  TAPE  3,30,  NBRAK,NDUP,NLNK,AN«W,ASMW,TIMF 
30  FORMAT  I1X,3I10,3F20,5) 

25  CONTINUE 
CALL  EXIT 
fcNO 
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APPENDIX  II 
COMPUTER  PROGRAM 
SIMPLE  CASE.  TYPE  II 
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TESSLER  ladder  polymer  degradation  fused  6  MEMBEREO  RING 
XEQ 

DIMENSION  Lfl00,200» 

FREQUENCY  7 » 1 *0,2 >  *8 (9  *  1 »0 J ♦ 23 ( 1 » 1 » 0 ) ,40 < 0* 1 #9 »  *42  1 0,1 *9 J 
1*4110*1,91 *43 10,1,9 1,4 1  2,0*1 » ,10 {9,1,0 J,21( 1,1,0 J,47C 0,1,91 
2,49 ( 0,1,9 > ,48  to, 1,9) ,50(0,1 ,9 » ,14(0,1 ,9* 

NBRAK=0 

N0UP=0 

NLNK=0 

DO  25  1X=1,30 
DO  24  1=1,200 
R=RAND1(Y >*19600. 

MR=R 

ML=NR/196+1 

NB=196*ML-NR 

7  lF|NB-78)8,8,4 

8  IF(L(ML,NB)-1)28,11,28 
28  L(ML,NB)=1 

NA=NB/2 

23  1F(2*NA-NB)40,41,40 

40  IFCl-L(ML,NB+78) )42,5,42 

42  lF(l-L(ML*NB+79) )24,5,24 

41  lF(l-L(ML,NB+78) )43,5,43 

43  lF(l-L(ML,NB+77) )24,5,24 

4  rF(NB-156)lO,10,14 

10  1F(L(ML,NB)-1)17,11,17 
17  L(ML,NB)=1 

NCsNB/Z 

21  IF (2*NC-NB)47, 48,47 

47  IF(l-L(ML,NB-78) )49,5,49 

49  IF(l-L<ML,NB-77) )24,5,24 

48  IF(l-L{ML,NB-78) )50,5,50 

50  IF(l-L(ML,NB-79) )24,5,24 
14  IH1-L(ML,NB)»  12,11,12 
12  L(ML,NB)=1 

NLNIC=NLNK+1 
GO  TO  24 

5  NBRAK»N8RAK+1 
GO  TO  24 

11  N0UP=NDUP+1 

24  CONTINUE 
BRK=NBRAK 

AMMW=100,/<100.+BRK) 

XI  =  IX 
0UPN=NDUP 

T1ME=L0GF (19600./ (19600, -200, *X  I +DUPN) ) 


ASMW=100./(100,+200,*Xl-0UPNj 

WRITE  OUTPUT  TAPE  3 , 30 ,NBRAIC  ,NDUP ,NLNK  ,  ANMW ,ASMW,T IME 
30  FORMAT(1X,3I10,3F20,5) 

25  CONTINUE 
CALL  EXIT 
END 

J1  ex-  -  IG7 

*J  -9(NZ5GA  7  -2  -E  -6  -  35-  -  97  95  --  94  78  ~  -5  -5  -7  —  4E  - 
•07IPCST  9  -G  5484  -1  -9  -5  -  2(X)  -  9  - 


070SQ000 

-070SOU01 

070_sflU02 
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TESSLER  LADDER  POLYMER  DEGRADATION  WHERE  BRAK  LNKS  CAM  YIELD  BRAK 
XEO 

DIMENSION  L(100,200l 

FREQUENCY  7J 1 .0*2 J *8 (9 » 1 *0 ) 1 9 f 0 . 1 ,9 ) *31 ( 1 *1 *81 1 0*1 ,9> *84(0*9.1 ) 
1 *93(0*1.9) *54(0 *1.9) ,69(9.1 .1) .42(0,1*9) .44 (0*9* 11*46(0*1 *9 > 

2. 48  (0*1 .9)  ♦41 1,0,1)  .32(9,1,0)  .6  (0.1,9 )  ,97  (1.1  ♦$()  *35(9.1.0) 

3*36 (0*9.1 ) *60(0,1 *9 ),62( 0,1 *9) *68 (9, 1,1) .20 (0,1, 9). 21 (0,9,1) 

4,91 (0*1,9), 27(0, 1,9) ,13(9,1*0) ,190 ( 1,1 ,9) ,191 (9,1, 0),107( 9, 1,0) 
5,109(9*1,0) ,  111(9,1*1 ) ,194 (9,1,0), IIB (9, 1,0) 

6,120(9*1*0) 

NBRAK=0 

N0UP=0 

NLNK=0 

DO  25  IX=1,30 
DO  24  1*1*200 
k=RAND1(Y)»19900. 

NR*R 

MLaNR/199+l 

NB*199*ML-NR 

7  IF(N8-66)8,8,4 

8  rF(L(ML,NB)^l)28, 11*28 

28  L{ML,NB)=1 

9  IF(l-L(ML,NB+66) )80,5,80 
C  OPPOSITE  BOND  CHECKED 

80  DO  51  rc:=l,6T 
NW*NB+133-IC 

31  1F(NW-133)5,81,81 

81  1F(1-L(ML*NW) )84,51,84 

51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

84  IF|IC-1)85*75,85 

85  10*1C-1 

DO  52  KA*l,IO 
NX=N3+66-KA 

93  IF(1-L(ML,NX) )53,5*53 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 
33  NYsNB-KA 

54  1F(1-LCHL,NY))52,5*52 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

52  CONTINUE 

75  DO  17  1C»1,67 

NW=NB-l-132+IC 
69  IF(NW-199)42*42*5 

42  1F(1~L(ML,NW))44,17,44 

17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 
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44  rF<IC-l)45*24»45 

45  10=tC-l 

DO  19  KA»1*I0 
NX=NB+66+ICA 

46  IF(1-L(ML*NX) »47»5»47 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

47  NYaNB+KA 

48  IF(1-L(ML*NV»}19»5»19 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

19  CONTINUE 
60  TO  24 

4  IF(NB-132)32»32»13 

32  IF(L(ML*NB)-1»33»11»33 

33  L(ML«NB)al 

6  IF(l-L(ML»NB-66n34,5»34 

C  OPPuaiTE  BOND  CHECKED 

34  DO  90  IC:=*1*67 
NW=NB+67-IC 

97  IF(NW-133)5*35*35 

35  If (L(ML*NW)-1)36.90,36 

90  continue 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  IF(IC-1137»38»37 

37  I0«IC-1 

DO  39  <A»1»I0 
NX*NB-66-KA 

60  IF(l-L(ML*NX))61»5*6i 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

61  NY=NB-KA 

62  IF(1-L(ML*NY»139»5»39 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUE 

38  DO  63  IC>=1»67 
NW*NB+66+IC 

68  iF(NW-199>20»20,5 

20  IF(1-L(ML*NWM21»63*21 

63  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

21  IF( IC-1 >22»24*22 

22  I0=IC-1 

DO  23  KA*l»tO 
NX»NB-66+KA 

91  1F( 1-L(ML*NX))26»5»26 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

26  NY=N8+KA 

27  IF(1-L(ML*NYM23»5»23 
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C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

23  CONTINUE 
60  TO  24 

13  IF(HML#NBJ-1)29»11#29 

29  LfML.NBI=l 
NLHK=NLNK+1 

DO  100  IC=1»67 

190  IF{NB-IC-133>104*102.102 
102  IP=NB-IC 

191  tFJL(ML»tPl-l)l05.100,105 
100  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  DO  106  1M=1»IC 
NX*NB-1W-132 

107  lF<LrML*NXT-l>108.104»104 

108  NY=NB-1N-66 

109  IFCL(ML*NY)-1J106*104*104 

106  CONTINUE 
GO  TO  24 

C  CHECKED  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  110  IC=1.67 

111  1F(NB+IC-199)112»112»5 

112  IP=NB+IC 

194  IF(L(ML.IP)--1)116.110.116 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

116  DO  117  IW=1»IC 
NX=NB+IW-r33 

118  IF{LCML»NX)-1J119,5,119 

119  NY*NB+IW-67 

120  IF(L{ML»NY)-1)117.5»117 

117  CONTINUE 

C  CHECKED  ALL  SIDELINKS  TO  THE  RIGHT 
GO  TO  24  ' 

5  NBRAK=NBRAK+1 
GO  TO  24 
11  NDUP=NDUP+1 

24  CONTINUE 
BRK=NBRAK 

ANMW=100./(100.+BRK) 

XIalX 

DUPNaNOUP 

T IMEaLOGF 1 19900,/ (19900. -200. »X I+DUPN ) } 
ASMW=100./(100,+200,*XI-DUPN) 

WRITE  OUTPUT  TAPE  3 .30 ,NBRAK .NDUP .NLNK ♦ANNW *ASMW»T IME 

30  F0RMAT(1X»3I10»3F20,5) 


25  CONTINUE 

call  exit 

END 
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APPENDIX  IV 
COMPUTER  PROGRAM 
COMPLEX  CASE.  TYPE  II 
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C8  TYPE  II  WHERE  CROSSLINK  BREAK  LEADS  TO  ADDITIONAL  SIDE  LINK  BREAKS 
*  XEO 

DIMENSION  L(100,200» 

FREQUENCY  7 ( 1 ,0,2 ) , 8 { 9 , 1 .O) *92 < 1 » 1 , 0 ) 0* 1 *9 ) *95 <0,1 *9» 

1,94 <0,1, 9) ,96(0,1,9) ,31 (1,1,9) ,8 1 ( 0 , 1 ,9 ) ,84 (0,9 , 1 ) ,73 ( 0, 1 ,9) 
2,54(0,1 ,9), 69 (9, 1,1) ,42(0,1 ,9) ,44 (0,9,1), 46 (Oil ,9), 48 (0,1, 9) 
3,4(1,0,1) ,32(9,1,0) ,100 ( 1,1 ,0) ,98(0,1,9) *101(0,1,9) ,99(0,1,9) 

4, 102 (0,1,9) ,97 (1,1,9) ,35(9,1,0) ,36(0,9,1) ,60(0,1,91,62(0,1,9) 
5,68 (9,1, 1 ),20{ 0,1,9) ,21(0,9,1) ,91 (0,1,9 ),27 (0,1,91,13 (9,1,0) 

6, 190 (1,1,9), 191 (9,1,0) ,107(9,1,0) ,109(9,1,6) 

7, 111(9, 1,1), 194 (9, 1,0), 118(9, 1,0), 120 (9, 1,0) 

NBRAK*0 

NDUP=0 

NLNK=0 

DO  25  IX=1,30 
00  24  1=1,200 
R=RAND1(Y)*19600, 

NR=R 

ML=NR/196+1 

NB=196*ML-NR 

7  IF(NB-78)8,8,4 

8  rF(L(ML,NB)-l)28*ll,28 
28  L(ML*NB)=1 

NA=NB/2 

92  1F(2*NA-NB)93,94,93 

93  IF(l-L(l5L,NB+78))95,5,95 

95  IF(l-L(ML,NB+79))80,5,80 

94  1F(1-L(ML,NB+78) 196,5,96 

96  IF(l-L(ML,NB+77))80,5,80 
C  BOTH  OPPOSITE  BONDS  CHECKED 

80  DO  51  lC=i,41 
NW=NB“NB/2+157-lC 

31  1F(NW-157)5,81,81 

81  1F(1-L(ML,NW) 184,51,84 

51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

84  IF(IC-1)85,75,85 

85  IO=(IC-l)*2 
DO  52  KA=l,IO 
NX=((NB+l)/2)*2+77-KA 

73  1F(1-L(ML,NX))53,5,53 
C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

53  NY=NX-78 

54  IF(1-L(ML*NY))52,5,52 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

52  CONTINUE 
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75  DO  17  IC=1*41 

NW=NB-NB/2+156+IC 
69  IF(NW-196 U2.42.5 
42  IF(1-L(ML»NW) )44»17*44 
17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

44  IF( IC-1 )45»24.46 

45  IO=(IC-l)*2 
DO  19  KA=l.IO 

NX  =  nNB+l  )/2)*2+78+IO 

46  IF(1-L(ML.NX) )47»5*47 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

47  NY=NX-78 

48  IF(1-L(ML*NY) )19.5»19 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
19  CONTINUE 
GO  TO  24 

4  IF(NB-156)32*32.13. 

32  IF(L{ML.NBr-l)33»11.33 

33  L(MLtNB)*=l 
NC=NB/2 

100  IF(2*NC-NB)98,99f98 

98  IF(l-L<ML»NB-78) )101 .5*101 

101  IF(l-L(ML.NB-77) 134,5.34 

99  I F ( 1 -L  <  ML .NB-78 ) ) 1 02 ,5,102 

102  IF{l-L(ML,NB-79) )34,5,34 
C  OPPOSITE  BONDS  CHECKED 

34  DO  90  IC=1*41 
NW=NB-NB/2+118-IC 

97  IF(NW-157)5,35*35 

35  IF(L(ML,NW)-1)36,90,36 
90  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  IF( IC-1 )37, 38,37 

37  10= ( IC-11»2 
DO  39  KA=l,IO 

NX=( (NB+1 )/2)*2-79-KA 

60  IF(1-L(ML,NX)161,5,61 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

61  inY=NX+78 

62  IF(1-L(ML,NY) )39»5,39 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUE 

38  DO  63  IC=1,41 
NW=NB-NB/2+117+TC 

68  IF (NW-196 )20»20,5 
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20  1F<1“L(ML»NW))21.63.21 
63  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

21  IFIIC-l )22, 24*22 

22  IO={IC-l)*2 
DO  23  KA=l.IO 

NX=( (NB+1 )/2)*2-78+KA 
91  IF{1-L<ML»NXJ)26.5*26 
C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

26  NY=NX+78 

27  IF(1-L(ML.NY> )23.5*23 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

23  CONtiNUE 
GO  TO  24 

13  1F(L<ML.NB)-1)29.11.29 
29  L(ML*NB)=l 
NLNK=NLNK+1 
DO  200  IC=1.41 

190  IF(NB-IC-157)104*202.202 
202  IP*NB-IC 

191  IF(L(ML»IP)-1)105.200.105 
200  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  IBA=2*IC 

DO  106  IW=1*1BA 
NX=2*(NB-156)-1-IW 

107  IF(L(ML»NX)-1)108. 104*104 

108  NY=NX+78 

109  IF(L(ML*NY)-1)106*104*104 

106  CONTINUE 
GO  TO  24 

C  CHECKED  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  110  IG=1*41 

111  IF<NB+IC-196)il2*112*5 

112  IP=NB+IC 

194  IF{L{HL*IP)-l>116*nO*116 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

116  IBA=2»IC 

DO  117  IW=1*IBA 
NX=2*(NB-157)+78+IW 

118  1F(LJML*NX)-1)119*5*119 

119  NY=NX-78 

120  iF{LIML*NY)-l)117*5*117 

117  CONTINUE 

C  CHECKED  ALL  SIDELINKS  TO  THE  RIGHT 
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60  TO  24 
5  NBRAK*NBRAX+1 
GO  TO  24 
11  NDUP»NDUP+1 

24  CONTINUE 
BRK=NBRAK 

ANMW=100./( lOO.+BRK) 

XI  =  IX 
DUPN=NDUP 

T  I ME*LOGF(  19600./ (19600. -200. *XI+DUPNn 
ASMW=100./( 100.+200.»XI-DUPN) 

WRITE  OUTPUT  TAPE  3 . 30 ,NBRAK »NDUP »NLNK » ANMW » ASMW . T I  ME 
3U  m>KMAT( IX .31 10.3F20.5 ) 

25  CONTINUE 


CALL  EXIT 
END 


)1  (X-  - 

*1  -9(NZ5GA 

(G7 

7  -2 

-E  -G  - 

3  5- 

-  97  95  — 

94  78  — 

070SOU00 

_5  _5  -7  —  ZlF  -4  -070SOU01 

*U7(P($T 

9  -G 

5484  -1 

-9  -5  - 

2(X> 

-  9  - 

ri70<;Q'J0? 
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APPENDIX  V 
COMPUTER  PROGRAM 
SIMPLE  Case,  type  i 
SINGLE  BOND  EVERY  SIXTH  BOND 
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CIS  EVERY  SIXTH  BOND  IS  SINGLE  CHAIN  BOND  CALCULATE  NAMM  CHANGE 
C  CJISE  I  SIMPLE  PROBABILITY  FOR  ALL  BONOS  BREAKING  IS  EQUAL 
DIMENSION  LC100,200J 

FREQUENCY  7 C 1 »0 , 2 »  *8 (9 . 1 #0 ) ♦ 91 0 ♦ 1 *9  I *41 1 »0 ♦ 1 J  » 1 0 { 9»1 *0 ) , 1 T (0 , 1 *9 J 
l*14I9tl.0> 

NBRAK=0 

NOUP=0 

NLNK=0 

00  810  MA=1»100 
00  820  IB=6»66*6 
820  LCMA»IB)=1 
810  CONTINUE 

00  25  IX=1»30 
DO  24  1=1*200 
R=RANDHY»*19900, 

NR=R 

ML=NR/199  +1 
NB=199*ML  -  NR 
7  IFINB-6618,8,4 

8  1F(L(HL*NBJ-1)28»11*28 
28  L(ML*NBJ=1 

9  IFa-L(ML*NB+66)  J24,  S*24 

4  IF(NB-132)10, 10*14 

10  IFa{ML»NB)-l)17»ll*17 

17  L(ML.NBI=1 

13  IF(l-HML»NB-66)  124*5*24 

14  IF(LrML*NB)-l)12»ll»12 

12  LtML»NBJ=l 

NLNK=NLNK+1 
GO  TO  24 

5  NBRAK=NBRAK  +1 
GO  TO  24 

11  NDUP=NDUP+1 

24  CONTINUE 
BRK=NBRAK 

ANMW=100./n00,+BRK» 

X1=1X 

DUPN=NOUP 

T 1ME=L0GF <18800. /( 18800, -200, *XI+OUPN » } 
ASMW=100,/<100,+200,*XI-0UPN) 

WRITE  OUTPUT  TAPE  3*30,  NBRAK »NDUP*NLNK*ANMW.ASMW»T1ME 
30  FORMAT  C1X,3I 10,3F20,5 ) 

25  CONTINUE 
CALL  EXIT 
END 
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APPENDIX  VI 
COMPUTER  PROGRAM 
SIMPLE  CASE.  TYPE  I 
SINGLE  BOND  EVERY  TWELFTH  BOND 
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C  CASE  1  SIMPLE-PROBABILITY  FOR  ALL  BONOS  BREAKING  IS  EQUAL 
C16  EVERY  TWELVTH  BOND  IS  SINGLE  CHAIN  BOND  -CALCULATE  NAMW  CHANGE 
DIMENSION  L( 100,200 J 

FREQUENCY  7(l,0,2j»8{9»l,0>,9(0,l,9)»4(l»0,l),l0<9»l»0»,13l0»l,9» 
1,1«I9»1»0> 

NBRAK=0 

NDUP=0 

NLNK=0 

DO  810  MA=1»100 
DO  820  1B=12»66»12 
820  L{MA,IB)=1 
810  CONTINUE 

DO  25  IX=1»30 
DO  24  1*1,200 
R=RANDltY >*19900. 

NR=R 

ML=NR/199  +1 
NB=199*ML  -  NR 
7  rF(NB-66)8»8,4 

8  1FIL(ML,NB)-1)28»11,28 
28  L<ML,NB)*1 

9  IF(1-L(ML»NB+66»)24,  5.24 

4  IF(NB-132 >10,10.14 

10  IF(L(ML.NB)-1>17.11,17 
17  L(ML.NB>=1 

13  IF(l-L(ML.N8-66))24,5.24 

14  IF(L{ML.N8)-1)I2.11.12 
12  L{ML,NB>=1 

NLNK=NLNK+1 
GO  TO  24 

5  NBRAK=NBRAK  +1 
GO  TO  24 

11  NOO|i*NDUP+l 

24  CONTINUE 
BRK=NBRAK 

AlilMW*100,  /( lOO.+BRK  > 

XI*IX 

OUPN=NOUP 

T 1 ME*LOGF ( 19400 . / 1 1 9400 . -2  00 . *X I +OUPN  > ) 

ASMW=100./« 100.+200,*XI-DUFN) 

WRITE  OUTPUT  TAPE  3.30,  NBRAK.NDUP.NLNK.ANMW.ASMW.TIME 
30  FORMAT  I1X,3110,3F20.5 > 

25  CONTINUE 
CALL  EXIT 
END 
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COMPUTER  PROGRAM 
SIMPLE  CASE.  TYPE  II 
SINGLE  BOND  EVERY  TWELFTH  BOND 
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C19  CASE  II  SIMPLE  EVERY  TWELVTH  BONO  IS  SINGLE  CHAIN  BOND 
C  PROBABILITY  FOR  ALL  BONDS  BREAKING  IS  EQUAL-CALCULATE  NAMW  CHANGE 
*  XEQ 

DIMENSION  LI100,200) 

FREQUENCY  7 ( 1 .0,2 ) » 8 ( 9 » 1 *0 ) , 23 ( 1 ♦ 1 ♦ 0 ) *40 1 0 » 1 ,9 ) »42 ( 0»1 *9» 

1,41 <0,1,9  1, 4310,1,91 ,4 (2.0.1 ) .10(9,1,0) ,21 (1.1,0 ),47(0.1,9I 
2, 49 (0,1, 9), 48 (0,1, 9), 50(0, 1,9), 14(0, 1,9) 

NBRAK=0 

NDUP=0 

NLNK=0 

DO  810  MA=1,100 
DO  820  18=12,78,12 
820  L(MA,IB)=1 
810  CONTINUE 

DO  25  IX=1,30 
Du  24  1=1 ,200 
R=RAND1(Y)*19600, 

NR=R 

ML=NR/196+1 

NB=196*ML-NR 

7  IF(NB-78)8,8,4 

8  IFtL(ML,NB)-l)28,ll,28 
28  L(ML,NB)=1 

NA=NB/2 

23  IF(2»NA-NB)40,41,40 

40  IF(l-L(ML,NB+78))42,5,42 

42  IF(l-L(ML,NB+79))24,5,24 

41  IF(1-L(ML,NB+78) 143,5,43 

43  IF( l-L(ML,NB+77) )24,5,24 

4  IF(NB-156)10,10,14 

10  IF(L(ML,N8)-1)17,11,17 
17  L(ML,NB)=1 

NC=NB/2 

21  IF(2*NC-NB)47,48,47 

47  IF(l-L(ML,NB-78))49,5,49 

49  IF(l-L(ML,NB-77) )24,5,24 

48  IF(l-L(ML,NB-78))50,5,50 

50  IF(l-L(ML,NB-79) )24,5,24 
14  1F(1-L(ML,NB) )12,11,12 
12  L»mL,NB)=1 

NLNK=NLNK+1 
GO  TO  24 

5  NBRAK=NBRAK+1 
60  TO  24 

11  NDUP=Nr)UP+l 

24  CONTINUE 
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BRK=NBRAK 

anmw=ioo./(ioo.+brk:) 

XI  =  IX 
DUPN=N0UP 

TIME=LOGF( 19000,/( 19000. -200 .*X I+DUPN ) ) 

ASMW=100./( 100.+200.*xl-DUPN  » 

WRITE  OUTPUT  TAPE  3 , 30 , NBRAK tNDUP .NLNK » ANWW ♦ ASMW .T IME 
30  FORMAT(1X,3I10,3F20.5) 


25  CONTINUE 
CALL  EXIT 
END 


)1  (X-  - 

♦J  -9(NZ5GA 

(67 

7  -2 

-E  -6  - 

3  5- 

-  97  95  —  94  78  ~ 

070SOU00 

-5  -5  -7  —  4E  -4  -070SOU01 

*U7(P($T 

9  -G 

5484  -1 

-9  -5  - 

2(X)  -  9  - 

07050002 
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APPENDIX  VIII 
COMPUTER  PROGRAM 
SIMPLE  CASE.  TYPE  II 
SINGLE  BOND  EVERY  TWENTY- FOURTH  BOND 
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XEQ 

DIMENSION  L(100»200)  4n{0.1*9l ♦4?(0#l»9l 

ia!ro!u,,lu:S:iUr;4:^;S ;?  io:,;i.;..nn.i^ 

2*49t0.1*9)»48r0,l*9) .50(G,1»9) .14(0*1»95 

MBRAK-0 

NDUP*0 

NLHX»0 

DO  810  NA=1»100 
00  820  1B=24*78»24 
820  HMAt!Bt  =  l 
810  CONTINUE 

DO  25  IX=1»30 


DO  24  1=1 *200 
R=RAND1{Y)»19600. 

NR*R 

ML *NR/ 196+1 
NB*196*ML-NR 

7  IF(NB-78)8»8.& 

8  IF<HMU»NB>-1)28*11»28 
28  L(ML»NB)*1 

NA=NB/2 

23  ih2«NA-NB)40,41.40 

40  1F(1-L(ML»NB+78)142»5»42 

42  IF{ 1-LIML»NB+79) )24.5«24 

41  IF(1-L(ML.NB+78) »43»5.43 

43  1F(1-L(ML»NB+77))24»5*24 
4  1F|NB-156)10.10»14 

10  !F<L{ML»NBI-1)17»11*17 
17  HML»NB»=1 
NC=NB/2 

21  IF(2*NC-NB)47»48*47 

47  IF(1-L<ML»NB-78))49,5»49 

49  IF(l-L(ML.NB-77) I24.5.24 

48  lFa-LtMLfNB-78)  )50t5»50 

50  1FI1-LIML*NB-79) I24»5»24 
14  lFfl-LIMt»NBH12»Tl*12 
12  HML»NBI»1 

NLMK=NLNK+1 


i-rt 


5  MBRAK»N8RAK+1 
GO  TO  24 
11  M0UP*N00P+1 
24  CONTINUE 
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BR<=NBRAK 

ANMW=100./(100.+8RK) 

XI  =  IX 
DUPN=NDUP 

TIME=LOGF{ 19300,/ t 19300.-200 .*X I+DUPN ) ) 
ASMW=100,/(100,+200.»XI-OUPN> 

WRITE  OUIHUT  TAPE  3,30,NPRAK*NDUP,NLNK»ANMW*A5MW*TIME 
30  FORMAT(1X»3I10,3F20.5) 

25  CONTINUE 
CALL  EXIT 
END 

)1  (X-  -  (67 

*J  -9(NZ56A  7  -2  -E  -G  -  35-  -  97  95  —  94  78  - 5  -5  -7  —  4E  - 

♦U7(P($T  9  -G  5484  -1  -9  -5  -  2(X)  -  9  - 


070se)tJ00 

-07OSOU01 

070SOU02 
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APPENDIX  IX 
COMPUTER  PROGRAM 
COMPLEX  CASE.  TYPE  I 
SINGLE  BOND  EVERY  SIXTH  BOND 
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C17  CASE  I  COMPLEX-EVERY  SIXTH  BOND  IS  SINGLE  CHAIN  BOND 
C  PROBABILITY  FOR  ALL  BONDS  BREAKING  IS  EQUAL  -  CALCULATE  NAMW  CHANGE 
*  XEO 

DIMENSION  LI100,200) 

FREQUENCY  7 1 1 *0 ,2 ) » 8 ( 9 . 1 ,0 ) . 9 ( 0 » 1 .9 ) ,31 ( 1  *  1 *9 ) * 81 C 0 ♦! *9 ) *841 0*9,I > 
1,9310,1,9  » ,54(0,1,9) ,69(9,1 ,1 > ,42(0,1,9) ,4410,9,1) ,46(0,1,9) 

2,48 (0,1,9), 4(1, 0,1) ,32(9,1,0) ,6 (0,1,9) ,97(1,1,9) ,35(9,1,0) 

3, 36  CO, 9, 1),60 (0,1, 9  I ,62(0,1,9) ,68 (9, 1,1), 20 (0,1,91,21 (0,9,1) 
4,91(0,1,9) ,27(0,1,9) ,13(9,1,0) ,190 ( 1,1,9) ,191 (9,1,0 ),107( 9,1,0) 
5,109(9,1,0) ,  111 (9,1,1 ) ,194(9,1,0) ,118(9,1,0) 

6,120(9,1,0) 

NBRAK=0 

NDUP*0 

NLNX=0 

DO  810  MA=1,100 
00  820  18=12,66,12 
820  L(MA,IB)=1 
810  CONTINUE 

DO  25  IX=1,30 
DO  24  1=1,200 
R=RAN01(Y)*19900. 

NR=R 

ML=NR/199+1 

NB»199*ML-NR 

7  IF(NB-66)8,8,4 

8  1F(L(ML,NB)-1)28,11,28 

28  L(ML,NB)=1 

9  1F( l-L(ML,NB+66) )80,5,80 
C  OPPOSITE  BOND  CHECKED 

50  DO  51  IC=1,67 
NW=NB+133-IC 

31  1F(NW-133)5,81,81 

81  1F(1-L(ML,NW) )84,51,84 

51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECI^EO 

84  1F( 1C-1)85,75,85 

85  10=1C-1 

DO  52  KA«l,IO 
NX=NB+66-KA 

93  1F|1-L(ML,NX))53,5,53 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

53  NY=NB-KA 

54  1F(1-L(ML,NY))52,5,52 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

52  CONTINUE 

75  DO  17  IC»1,67 


40 


ML  TDR  64-151 


NW=N8+132+IC 
69  IFtNW-199K2»42.5 

42  IF(1-L(ML»NW)J44*17,44 

17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

44  IF( IC-1 )45,24.45 

45  10=IC-1 

DO  19  KA=l.IO 
NX=NB+66+KA 

46  1K1-L(ML.NX))47»5.47 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

47  NY=NB+KA 

48  1F{ 1-L(ML»NY) )19»5»19 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

19  CONTINUE 
60  TO  24 

4  IF(NB-132)32.32«13 

32  IF(L(ML»NB)-1J33»11.33 

33  L(ML»NB>=1 

6  IF(1-L(ML.NB-66H34,5*34 

C  OPPOSITE  BONO  CHECKED 

34  DO  90  IC*1»67 
NW=NB+67-IC 

97  IF(NW-133)5.35.35 

35  IF(L(ML»NW)-1)36*90,36 
90  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  IFI  IC-1  )37, 38*37 

37  I0=IC-1 

DO  39  KA=1»I0 
NX=NB-66-KA 

60  IF(l-L(ML»NXn61*5»61 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

61  nt=nb-Ka 

62  IF (1-L(ML*NY») 39*5*39 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUt 

38  DO  63  IC=1*67 
NW=NB+66+IC 

6b  IF(NW-199)20*20*5 

20  IF(  l-L(ML*NWn21*63,21 

63  CONTINUE 

C  ALL  crosslinks  TO  THE  RIGHT  CHECKED 

21  1F( IC-1)22*24*22 

22  I0*IC-1 

DO  23  KA=1*I0 
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NX=NB-66+KA 

91  1FC1-HML»NX)  J26*5»26 
C  ALL  OPPOSITE  SIDE  LINKS  rHECKEO 
Zb  NT=nIJ+KA 

27  IF(l-LCML»NYn23»5»23 
C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

23  CONTINUE 
GO  TO  24 

13  IFIHML*NB)-1)29»11.29 

29  HML,NBJ  =  1 

NLNK=NLNK+1 
DO  100  IC=1,67 

190  IF(NB-1C-133)104»102.I02 
102  IP=NB-IC 

191  IF(L(ML»IP)-1)105,100»105 
100  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  DO  106  IW=1»IC 
NX=NB-IW-132 

107  IF(L<ML»NX>-1)108»104*104 

108  NY=NB-IW-66 

109  IFIL<ML»NY)-1)106.104.104 

106  CONTINUE 
GO  TO  24 

C  CHECKED  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  no  IC  =  1»67 

111  IF(NB+IC-I99)112»112»5 

112  IP=NB+IC 

194  IF(L(ML»IP)-1)116.110»116 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

116  DO  117  IW=1.IC 
NX*NB+IW-133 

118  IF(L(ML»NX)-1)119»5,119 

119  NY=NB+IW-67 

120  IF{L(ML*NY)-1)117»5,117 

117  CONTINUE 

C  CHECKED  ALL  SIDELINKS  TO  THE  RIGHT 
GO  TO  24 

5  NBRAK=NBRAK+1 
60  TO  24 
11  NDUP=NDUP+1 

24  CONTINUE 
BRK=NBRAK 

ANMW=100./( lOO.+BRK ) 

XI  =  IX 


DUPN=NDUP 

TIME=LOGF( 18800. /( 18800.-200.*XI+DUPN ) > 

ASMW=100,/( 100.+200.*Xr-DUPN  t 

WRITE  OUTPUT  TAPE  3 ,30 ,NBRAK .NDUP .NLNK »ANMW»ASMM,T IME 
30  FORMATI1X,3I10,3F20.5) 

25  CONTINUE 
CALL  EXIT 
END 
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APPENDIX  X 
COMPUTER  PROGRAM 
COMPLEX  CASE.  TYPE  I 
SINGLE  BOND  EVERY  TWELFTH  BOND 
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Cl6  CASE  I  COMPLEX  EVERY  TWELVTH  BOND  IS  SINGLE  CHAIN  BONO 

PROBABILITY  FOR  ALL  BONDS  BREAKING  IS  EQUAL  -  CALCULATE  NAMW  CHANGE 
XEO 

DIMENSION  L( 100,2001 

FREQUENCY  7 ( 1 ,0 ,2 ) ,8 { 9 , 1 ,0 ) , 9 ( 0 , 1 ,9  I  *31 f 1  *  1 ,9 ) » 81 < 0 ,1 ,9 ) , 84 1 0,9,1 1 
1,93 (0,1, 9 >,54(0,1,9) ,69 ( 9, 1  , 1 )  ,42 ( 0 , 1 ,9 ) ,44 ( 0*9 ♦ 1 ) ,46 ( 0 , 1 ,9 ) 
2,48(0,1,9) ,4 (1,0,1 ) ,32(9, 1,U) ,6(0, 1,9) ,97(1 ,1,9) ,35(9,1 ,0) 
3,36(0,9,1 ) ,60(0,1,9) ,62(0,1 ,9) ,68(9,1,1) ,20(0,1,9) ,21(0,9,1 ) 

4, 91 (0,1, 9 ),27(0, 1,9) ,13(9,1 ,0) ,190 ( 1,1,9) ,191(9,1,0) ,107(9,1 .rn 
5,109(9,1,0),  111 (9,1,1  )  ,194(9,1 ,0) ,118(9,1 ,0) 

6,120(9,1,0) 

NBRAK*0 

NDUP*0 

NLNK=0 

DO  810  MA=1,100 
DO  820  IB=6,66,6 
820  L(MA,IB)=1 

810  continue 

UO  25  IX=1,30 
DO  24  1=1,200 
R=RANDlfY)*19900. 

NR*R 

ML=NR/199+1 

NB»199*ML-NR 

7  rF(NB-66)8,8,4 

8  IF(L(ML,NB)-1)28,11,28 

28  L(ML,NB)=1 

9  IF(l-L(ML,NB+66) >80,5,80 
C  OPPOSITE  BOND  CHECKED 

80  DO  51  IC=1,67 
NW=NR+133-IC 

31  IF(MW-133)5,8l,8l 

81  IF(1-L(ML,NW) )84*5l,84 

51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

84  IF(1C-1)85,75,85 

85  IO=IC-l 

DO  52  KA=1,I0 
NX»NB+66-KA 

93  IF(1-L(I4L,NX))53,5,53 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

53  NY=NB-KA 

54  IF(1-L(W.*NY))52,5,52 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

52  CONTINUE 

75  DO  17  1C=1,67 
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NW=NB+132+rC 
69  1F(NW-199)A2»42»5 

42  IF(1-L«ML»NW) )44»17.44 
17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

44  IF( IC-1 )45»24»45 

45  10=1C-1 

DO  19  KA=l»IO 
NX=NB+66+KA 

46  IF{1-L(ML»NX) )47*5*47 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

47  NY=NB+KA 

48  rF(i-L(ML.NY) )19»5»i9 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

19  CONTINUE 
GO  ro  24 

4  IF(NB-132 )32*32»13 

32  IF(L(«L*NB)-1)33.11,33 

33  L(ML*N8)=1 

6  IF{l-L(ML.NB-66) )34»5»34 

C  OPPOSITE  BOND  CHECKED 

34  DO  90  IC=1*67 
NW=NB+67-IC 

97  IF (NW- 133) 5 ♦35.35 

35  IF{L(ML.NW)-1)36.90,36 
90  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  IF( IC-1 )37.38.37 

37  IO*IC-l 

DO  39  KA=l.IO 
NX=NB-66-KA 

60  IF(1-L(ML.NX) )61.5.61 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

61  NY=NB-KA 

62  IF(1-L(ML.NY))39.5.39 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUE 

38  DO  63  IC»1.67 
NW=N8+66+IC 

68  IF(NW-199)20.20,5 

20  IF{1-L(ML.NW))21.63,21 

63  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

21  IF(IC-1)22. 24.22 

2?  I0=IC-1 

DO  23  KA»1.I0 
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NX=NB-66+KA 

91  IFC1-L(ML*NX))26»5.26 

C  ALL  OPPOSITE  SIDE  LINKS  CTfECKEO 

26  NY=NB+KA 

27  IF(1-L(ML»NY)I23»5»23 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
23  CONTINUE 
GO  TO  24 

13  IFILJMLfNBI-1 J29.ll.29 
29  LIML.NBI=1 
NLNK=MLNK+1 
DO  100  IC*1.67 

190  IFfNB-IC-133J 104.102.102 
102  1P=NB-IC 

191  IF(L(ML»IPJ-1I105. 100.105 
100  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  DO  106  IW=1.IC 
NX=NB-1W-132 

107  IF(LCML»NXJ-1)108.104.104 

108  NY=NB-IW-66 

109  IF|L(ML.NY)-1>106.104.104 

106  CONTINUE 
GO  TO  24 

C  CHECKED  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  110  IC=1.67 
ni  IF(NB+IC-199)112»112.5 
112  IP=NB+IC 

194  IF(L(ML.IPJ-l>116.il0.116 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

116  DO  117  IW*1.1C 
NX=NB+IW-133 

118  IF(L(ML.NX}-1)119»5,119 

119  NY=NB+IW-67 

120  IFrL(ML.NYJ-l»117.5.117 

117  continue 

C  CHECKED  ALL  SIDELINKS  TO  THE  RIGHT 
GO  TO  24 

5  NBRAK=NBRAK+1 
GO  TO  24 
11  NDUP=NDUP+1 
29  CONTINUE 
BRK=NBRAK 

ANMW=100,/I100.+BRKI 
XI  =  1X 


DUPN=NDUP 

TIME»LOGF 119400./ (19400. -200. *X1+DUPN) J 
ASMW=100./( 100,+200.*XI-0UPNJ 

write  output  tape  3.30,NBRAK.NDUP.NLNK.ANMM.ASMW.TIME 
30  FORMAT(1X.3I10.3F20.5J 

25  CONTINUE 
CALL  EXIT 
END 
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APPENDIX  XI 
COMPUTER  PROGRAM 
COMPLEX  CASE.  TYPE  II 
SINGLE  BOND  EVERY  TWELFTH  BOND 
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C21  CASE  II  complex  EVERY  TWELVTH  BOND  IS  SINGLE  CHAIN  BOND 

PROBABILITY  FOR  ALL  BONUS  BREAKING  IS  EQUAL  -  CALCULATE  NAMW  CHANGE 
XEQ 

DIMENSION  LC100,200> 

FREQUENCY  7Cl»0,2J.8(9.1»0)»92(ltl»0),93(0,l»9>,95(0,l,9) 
lt94l0,l»9)»96C0,l,9>,31Cl,1.9)»81(0,l,9>»84(0,9,l),73(0,l,9J 
2»54(0»1,9) »69C9»1»1»  *42(0,1,9) ,44 (0,9,1), 46 (0,1,9), 48 (0,1, 9) 

3,4 (1,0,1) ,32(9,1,0) ,100 ( 1,1,0), 98 (0,1,9) ,101(0,1,9) ,99(0,1,9) 
4,102(0,1,9) ,97(1,1*9) ,35(9,1,0) ,36(0,9,1)  ,60 ( 0 , 1  ,9  )  ,62 ( 0 ,1  ,9 ) 
5,68(9,1,1 ) ,20 (0,1, 9), 21 (0,9,1),  91 (0,1, 9), 27 (0,1, 9), 13(9, 1,0) 
b, 190 (1,1, 9), 191 (9, 1,0), 107 (9, 1,0), 109 (9, 1,0) 

7, 111(9, 1,1), 194 (9, 1,0), 118 (9, 1,0), 120 (9, 1,0) 

NBRAK=0 

NDUP=0 

NLNK=0 

DO  810  MA=1,100 
DO  820  18=12,78,12 
820  L(MA,IB)=1 
810  CONTINUE 

DO  25  IX=1,30 
DO  24  1=1,200 
R=RAND1 (Y )*19600, 

'  NR=R 

ML=NR/196+1 

NB=196*ML-NR 

7  IF(N8-78)8,8,4 

8  IF(L(ML,NB}-1)28,11,28 
28  L(ML,NB)=1 

NA=NB/2 

92  IF(2*NA-NB)93,94,93 

93  IF(l-L(ML,NB+78) )95,5,95 

95  IK  l-L(ML*NB+79)  )'80,5,80 

94  IF(l-LIML*NB+78) )96,5,96 

96  IF( l-L(ML,NB+77) )80,5,80 
C  BOTH  OPPOSITE  BONDS  CHECKED 

80  DO  51  IC=1,41 
NW=NB-NB/2+157-IC 

31  IF(NW-157)5,81,81 

81  IF(1-LIML,NW) >84,51,84 
51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

84  IF(  IC-1 )85,75,85 

85  IO=(IC-l)*2 
DO  52  KA=l,IO 

NX={ (NB+1 )/2)*2+77-KA 
73  IF(1-L(ML,NX))53,5,53 
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C  ALL  OPPOSITE  SIDE  LINKS  CHECKFO 

53  NY=NX-78 

54  TF(1-L(ML*NY) )52»5.52 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
52  CONTINUE 
75  DO  17  IC=1»41 

NW=NB-NB/2+156+IC 
69  IF {NW-196  )42,42»5 
42  IF(1-L(ML»NW))44»17,44 
17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  R I GH I  CHECKED 

44  IF(IC-1 )45»24,45 

45  IO=(IC-l)*2 
DO  19  KA=l.IO 

NX=( (NB+1 )/2)*2+78+IO 

46  IF(1-L(ML»NX) )47»5»47 

C  all  OPPOSITE  SIDE  LINKS  CHECKED 

47  NY=NX-78 

48  IF(  1-L(ML  »NY')  )  19»5.i  Q 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
19  CONTINUE 
GO  TO  24 

4  IF(NB-156)32,32,13 

32  IF(L(ML»NB)-1 )33.11»33 

33  L(ML,NB)=1 
NC=NB/2 

100  IF(2*NC-NB)98*99*98 

98  IF(l-L(ML»NB-78) ) 101,5»101 

101  IF(l-L(ML»NB-77) >34*5,34 

99  IF(l-L(ML,rjB-78)  >102, 5*102 

102  IF( l-L(ML*NB-79) >34*5*34 
C  OPPOSITE  BONDS  CHECKED 

34  DO  90  IC=1*41 
NW=NB-MB/2+118-IC 

97  IF(NW-157>5*35*35 

35  IF(L(ML*NW>-1>36*90*36 
90  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  IF( IC-1 >37*38*37 

37  IO=(IC-l>*2 
DO  39  KA=l*IO 

NX= ( ( NB  +  1 > /2  >  *2-79-KA 

60  IF(1-L(ML*NX> >61,5*61 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

61  NY=NX+78 

62  IF(1-L(ML*NY> >39*5*39 
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C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUE 
38  DO  63  IC=1»41 

NW=NB-NB/2+117+IC 
68  IFINW-196120,20,5 

20  IF(l-HML»NWn21»63»21 
63  CONTINUE 

C  ALL  CROSSLINKS  lu  IHh  RIGHT  CHECKED 

21  IFC IC-1»22»24,22 

22  io=nc-ii*2 

DO  23  KA=1»I0 
NX=l  INB+l)/2 J*2-78+KA 
91  IF{l-L(ML,NXn26»5.26 
C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

26  NY=NX+78 

27  IFU-LCML»NYn23»5»23 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

23  CONTINUE 
GO  TO  24 

13  IFCL(ML»NBJ-1)29»11»29 
29  L(ML.NB)=1 
NLNK=NLNK+1 
DO  200  IC=1*41 

190  IF(NB-IC-157J104»202»202 
202  IP=NB-IC 

191  IF(L(ML.IP)-1)105.200,105 
200  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  IBA=2*1C 

DO  106  IW=1»IBA 
NX=2*f NB-156)-1-IW 

107  IF CL(ML»NX)-1I 108.104»104 

108  NY=NX+78 

109  IFILCMLjNY J-1) 106»104»104 

106  CONTINUE 
GO  TO  24 

C  checked  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  110  IC=1.41 

111  IF(NB+IC-196)112»112»5 

112  IP=NB+IC 

194  IF(LIML»IP»-l)ll6»110,iib 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 
116  IBA=2*IC 

DO  117  IW=1»IBA 
NX=2*INB-157 J+78+IW 
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118  IF(L(ML,NX)-1)119,5,119 

119  NY=NX-78 

120  IF(L(ML.NY)-1 )117»5.117 
117  CONTINUE 

C  CHECKED  ALL  SIDELINKS  TO  THE  RIGHT 
GO  TO  24 
5  NBRAK=NBRAK+1 
GO  TO  24 
11  NDUP=NDUP+1 
24  CONTINUE 
BRK=NBRAK 

ANMW=100./( lOO.+BRK) 

XI  =  IX 
DUPN=NDUP 

TIME=LOGF(19000./(19000.-200.*XI+DUPN) ) 

ASMW=lOO./( 100.+200.*XI-DUPN) 

WRirE‘'DUTPUT  TAPE  3 . 3U  ,  NBRAK  .NDUP  .NLNK.  ANMW  .  ASMW  *  T  I  ME 
30  FORMAT(1X.3I10»3F20.5) 


25  CONTINUE 
CALL  EXIT' 
END 
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C22  CASE  11  COMPLEX  -  EVERY  TWENTY  FOURTH  BOND  IS  SINGLE  CHAIN  BONO 
PROBABILITY  FOR  ALL  BONOS  BREAKING  IS  EQUAL  -  CALCULATE  NAMW  CHANGE 
XEO 

DIMENSION  LI  100*200  I 

FREQUENCY  71 1 .0,2 1  * 8 1 9 . 1 *0 >  *92 1 1 » 1 *0 ) ,93 1 0 ♦ 1 #9 1 *95 (0*1 *9) 

1, 9410,1*9 1,9610, 1,9 > ,311 1,1*9 1 *81 <0,1,9 ),8A(0,9,1) ,7310,1,9) 

2 *541 0*1 *9) *69 1 9, 1,1 1,42  1 0,1 *9) *44 <0,9*1), 46 <0,1 ,9) ,48(0,1,9) 
3*411,0*1) *32(9,1,0) *10011*1,0) ,9810,1*9) *101 1 0,1,9 ),99( 0,1,9) 

4 ,102  CO, 1*9) *97 1 1,1 *9) ♦ 35 ( 9 , 1 ,0 ) , 36 ( 0 ,9 *  1 ) ,60  1 0 , 1 ,9 ) ,62 1 0 , 1 ,9 ) 

5 *68 1 9* 1,1)  *20 10,1,9)  ,21(6,9,1) ,91(0,1,9) ,2710,1,9) ,1319,1,0) 
6,190(1,1,9), 191 (9*1,0) ,107(9,1,0) ,109(9,1,0) 

7*111 (9,1,1) *194(9, 1,0) *118(9,1*0) *120(9,1,0) 

NBRAK=0 

N0UP=0 

NLNK=0 

00  810  MA*1,100 
DO  820  16*24,78,24 
820  L(MA,1B)*1 

810  CONTINUE 

DO  25  IX*1,30 
DO  24  1*1,200 
R*RAND1 (Y) *19600. 

NR=R 

ML»NR/196+1 

N6*196*ML-NR 

7  IF(NB-78)8,8,4 

8  1F(L(ML,NB)-1)28*11,28 
28  L(ML*NB)=1 

MA*NB/2 

92  IF(2*NA>-NB)93,94*93 

93  lF(l-L(ML,NB+78))95,5,95 

95  lF(l-LlML,NB+79) )80,5»80 

94  lF(l-LIML,NB+78) )96*5,96 

96  IF(l-L(ML,N8+77))80*5,80 
C  BOTH  OPPOSITE  BONOS  CHECKED 

80  DO  51  IC=1*41 
NW=NB-NB/2+157-1C 

31  1F(NW-157)5,81,81 

81  IF(1-L(ML,NW))84,51,84 
51  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

84  IF(1C-1 )85,75,85 

85  I0=(1C-1)*2 
DO  52  KA=l,IO 
NX*l(NB+l)/2)#2+77-KA 

73  IFI1-L(ML,NX))53*5,53 
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C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

53  NYaNX-78 

54  IF(1-L(ML»NY) J52»5»52 

C  ALL  adjacent  side  LINKS  CHECKED 
52  CONTINUE 
75  DO  17  IC=1.41 

NW=NB»NB/2+156+IC 
69  IF(NW-196 (42»42»5 
42  IFfl-L(ML»NW))44»l7»44 
17  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

44  IFf IC-1 )45,24.45 

45  1U=(1C-1>*2 
DO  19  KA=l,ro 

NX={ (NB+1)/2)»2+78+I0 

46  IFn-L(ML»NX))47»5»47 

C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

47  NYaNX-78 

48  IF(1-L(ML»NY)  )19»5*19 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
19  CONTINUE 


GO  TO  24 

4  IF(NB-156l32.32»13 

32  1FIL(ML»NB)-1)33,11,33 

33  L(ML»N8>*1 
NCaNB/2 

100  IF(2*NC-NB)98,99»98 

98  lF(i-L(ML.NB-78) )101,5.101 

101  IF(1-L(ML»NB-77) )34,5»34 

99  IF( 1-L(ML»NB-78) ) 102»5*102 

102  IF (1-L( ML ♦N8~79)) 34,5,34 

OPPOSITE  BONDS  CHECKED 

„  34  DO  90  IC*1,41 

NW*NB-NB/2+118-IC 
97  1F(NW-157)5,35,35 
_35  IF(L(ML,NW)-1)36,90,36 
90  CONTINUE 

ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

36  1F(  iC-U37,38,37 

37  IO=.<IC-H#2 
DO  39  KA=1,10 

NXa(  (NB+1  J/2)»2--T9-KA 

60  IF(1-L(ML,NX»)61,5,61 

all  opposite  side  links  checked 

61  NY=NX+78 


62  IF(l-L(ML,NYn39»5,39 
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C  ALL  ADJACENT  SIDE  LINKS  CHECKED 
39  CONTINUE 
38  DO  63  IC=1»41 

NW=NB-NB/2+ll7+lC 
68  IFiNW-196 >20,20,5 

20  1F(1-L<ML»NW))21,63,21 
63  CCWTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 

21  IF<IC-1)22, 24,22 

22  IO=nC-lJ*2 
DO  23  KA=1,10 

NXa ( J  NB+1 » n >  *2-78+K A 
91  IF(1-HML»NXJ  I26»5»26 
C  ALL  OPPOSITE  SIDE  LINKS  CHECKED 

26  NY«NX+78 

27  IFri-LCML,NYH23,5,23 

C  ALL  ADJACENT  SIDE  LINKS  CHECKED 

23  CONTINUE 
GO  TO  24 

13  IFfL{ML#NBI-l>29»ll,29 
29  LtML#NB)»l 
NLNKoNLNK+1 
DO  200  rc*i»4i 

190  IF(N8-rC-157) 104,202,202 
202  IP»NB-IC 

191  IF(L(ML,IP>-1I105,200,105 
200  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  LEFT  CHECKED 

105  IBA=2*1C 

DO  106  IW«=1,1BA 
NX=2*iNB“156J-l-Hif 

107  IF(HML,NX)-1  >108,104,104 

108  NY«NX+78 

109  iFfLCML,NY)-l)106,l04,104 

106  CONTINUE 
GO  TO  2* 

C  CHECKED  ALL  SIDELINKS  TO  THE  LEFT 
104  DO  110  IG=1,41 

111  IF(NB+IC-196ni2,112,5 

112  IP»NB+IC 

194  IF|LIML,lPJ-mi6,110,H6 

110  CONTINUE 

C  ALL  CROSSLINKS  TO  THE  RIGHT  CHECKED 
116  I8A=2*IC 

DO  117  llif  =  l,lBA 
NX»2*INB-157>+78+lW 
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118  IF<L(ML*NX)-1)119*5»119 

119  NY=NX-78 

120  IF(L(ML»Ny)-l >117.5.117 
117  CONTINUE 

C  CHECICED  ALL  SIDELINKS  TO  THE  RIGHT 
GO  TO  24 
5  NBRAIC=NBRAK+1 
GO  TO  24 
11  NDUP=NDUP+1 
24  CONTINUE 
BRIC=NBRAK 

ANMW*100./( 100.+BRK> 

XI  =  IX 
DUPN=N0UP 

T I ME=L0GF{ 19300./ (19300, -200. *XI+0UPN) ) 

ASMW*100,/( 100,+200.»XI-DUPN) 

WRITE  OUTPUT  TAPE  3 . SO.NBRAK .NDUP.NLNK .ANMW.ASMW.T IME 
30  FORMAT(1X.3I10,3F20.5) 


25  CONTINUE 
CALL  EXIT 
END 
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